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Dependence of hopping conductance on temperature and voltage for an ensemble of modestly 
long one-dimensional wires is studied numerically using the shortest-path algorithm. In a wide 
range of parameters this dependence can be approximated by a power law rather than the usual 
stretched-exponential form. Relation to recent experiments and prior analytical theory is discussed. 

PACS numbers: 72.20.Ee, 73.63.Nm 
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Localization by disorder is a common cause of insu- 
lating behavior of low-dimensional electron systems. If 
the system size greatly exceeds the localization length, 
the transport at low voltages V and temperatures T is 
governed by the variable-range hopping (VRH) 1 . The 
electric current I{V, T) and the conductance G{V, T) = 
I{V,T)/V have a (stretched) exponential behavior. For 
example, in the Ohmic regime 

Go = G(0,T)cxexp[-(T,/T)^] , < /i < 1 • (1) 

Over the last decade, observations of different laws. 



(2) 
(3) 



have been reported in systems as diverse as carbon nan- 
otubes [H], InSb [TO] and GaAs [Hill] quantum wires, 
NiSes whiskers [13], polymer nanofibers [T?] [IS], inor- 
ganic [9] [16] and organic nanowires [17] , as well as poly- 
mer films [IHl [TO]- The coefficients a, /3, and 7 vary 
among different materials and different samples of the 
same material. 

A five-parameter formula frequently used to fit the ex- 
perimental data is 



/ = coT"+i sinh 
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For 7' = 7 the asymptotic behavior of /(V, T) is given by 
Eqs. (|2| and ([3|. Agreement with Eq. Q was advocated 
as evidence for tunneling into Luttinger liquid (LL) pU] 
— a one-dimensional (ID) system with nonperturbative 
interaction effects. (For strong interactions the LL can 
also be modeled as a ID Wigner crystal [21].) In this pic- 
ture, the system contains a tunneling barrier, e.g., a poor 
contact, but is otherwise clean and free of localization. 
The power-laws are due to renormalization of this bar- 
rier by many-body effects. However, there is a problem 
with this interpretation. The actual calculations [221(25] 
within the LL model give a — [3 and 7 = 7' = !, which 
is not always consistent with the parameters of the em- 
pirical fits (a notable exception is Ref . ITTj) . 



Another reason to doubt the relevance of the LL effects 
in some of these experiments is the fact that the systems 
studied are neither perfectly clean nor truly ID. They 
are, typically, collections of many parallel ID channels, 
whose total number ranges from several hundred to many 
thousands, each containing multiple impurities. 

In this Letter we show that in such quasi- ID systems 
the conventional mechanism of transport, which is the 
VRH, can also lead to Eqs. ([2])-(|4]). This is because at low 
enough T the hopping length is not much smaller than 
the length L of the wires. In this case, the VRH conduc- 
tance deviates from the usual formula, Eq. ([ij . The hop- 
ping is dominated by hopping paths that consist of a few 
approximately equidistant hops [2^1125] . Although rare, 
such configurations can always be found in a sample if the 
number of channels is large enough. Hence, despite meso- 
scopic fluctuations that accompany rare events, G{V,T) 
can be a smooth quasi-power-law function. 

Our results are illustrated in Fig. [l] They are obtained 
numerically following the approach used in our previous 
work [30] . with some improvements. In the inset of Fig. [l] 
we show a set of I-V curves computed for a set of fixed 
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FIG. 1. (Color online) Main panel: Collapse onto the "uni- 
versal curve" of Eq. (4) (solid line) obtained by plotting the 
numerical results for L = 30 and a = 4 (symbols) as j/y^+i 
vs. V/T, using a — 1.75, = 1.1, and 7 = 1. The tempera- 
tures are listed in the legend in the units of A = ATo- Inset: 
Same data plotted as I vs. V. 
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T. In the main panel, we collapse them onto a single 
"universal" curve described by Eq. Q. 

Let us compare the quality of the data collapse to those 
in the aforementioned experiments [51 riSHTOj . The range 
of T shown in Fig.[l]is a factor of 30. In the experiments, 
Eq. ([2]) rarely spans more than one decade in T. The 
range of V, where the non-Ohmic conductance follows 
the "universal" curve in the experiments, is usually less 
than a decade. In our case, it is wider than one decade. 
Still, the dependences that we find numerically are not 
true power-laws. If we look at wider ranges of V and T, 
the deviations are seen. Therefore, our numerical results 
for the VRH transport, just like the experiments, demon- 
strate only the apparent power-law behavior (APLB) re- 
stricted to a certain parameter range. 

In our calculations, this range is located near the in- 
flection point of the curve InGn vs. InT, see Fig. [2] 
Near the corresponding temperature Tinf the curve can 
be approximated by a straight line with a certain slope 
a, in agreement with Eq. ([2]). Further analysis, following 
Refs. which is discussed below, leads to analytical 

estimates 
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Here g is the density of states and q's are coefficients 
of the order of unity. In comparison, our simulations 
give a = 1.75, (3 — 1.1, and 7 = 1 for L/a — 7.5. For 
L/a = 12.5, we get a = 2.4, /3 = 1.7, and 7 = 0.6. This 
implies ci « 1.1, C2 ~ 0.85, C3 w 1.2, and C4 « 0.4. 

Our numerical results are comparable with typical ex- 
perimental numbers. They are also consistent with the 
observed trend that longer and more disordered wires 
produce larger a and /3 but smaller 7. A more detailed 
comparison would require taking into account particular- 
ities of a given set of samples beyond our generic model. 
Due to individual variations in the nature of disorder and 
the parameters of electron-phonon coupling, a and /3 may 
acquire additional corrections of the order of unity. 

Let us now give more details regarding the calcula- 
tions. We consider a system of localized states with 
random energies Si (Fig. [s]) distributed according to the 
Poisson distribution with the average energy separation 
A = 3. We treat all £i as constants, independent of the 
applied current. This is justified if electron interactions 
are weak [3D] . The x-spacing between the sites is taken to 
be unity, so that the density of states is g — 1/A = 1/3. 
The localization length a is chosen to be 4. To avoid 
repeated negative signs in the formulas, we take the elec- 
tron charge e = 1 to be positive. 

Let r]i be the electrochemical potential of site i, then 
the net current from site i to another site j is given by |31| 




FIG. 2. (Color online) The Ohmic conductance vs. temper- 
ature for L — 30 (upper curve) and L = 50 (lower curve) 
with a = 4. The dashed lines serve to illustrate the apparent 
linearity of the curves near their inflection points (dots). 
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where = F^ = Go |ei— £j|, and parameter Gq 

of dimension of conductance is related to electron-phonon 
coupling [TJ [28l [32] . The problem is to determine lij and 
rji that satisfy the current conservation. 

The analysis is simplified by the conventional assump- 
tion that the transport is dominated by a single optimal 
path of least resistance. Within this approximation, the 
current does not branch, i.e., lij — I in each link of the 
path. The total voltage drop V across the sample is the 
sum of voltage drops Ty^ — rjj on the links. One can de- 
termine the optimal path by finding the sequence of sites 
that gives the smallest V for a given I |30) . 

Let us define auxiliary variables ui = \n{TGo/I) and 



S = exp 

-L i 'i 
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Solving Eq. Q for rjj, we obtain 



= Tin 



(9) 



(10) 



Unlike the internal hops, transitions between the source 
electrode and the first site inside the sample (as well as 
the last site and the drain) do not require phonons. This 
can be accounted for by using 
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(11) 



in lieu of S. Here Gc is determined by the tunneling 
transparency of the contact between the sample and the 
electrode. We choose a representative value Gc = 4Go. 



Note that the numerator of Eq. (10 1 must be positive. 
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which sets a hmit on the maximum current that can flow 
between sites i and j. 

To find the optimal path, we use a modified Dijkstra 
algorithm [30] , in which the "cost" of reaching site j start- 
ing from the source equals —r]j. The latter is calculated 
recurrently using Eq. ([s]). For each disorder realization 
the resistance V/I is random and by running the simula- 
tions many times we can compute its probability distri- 
bution. Taking the average over the latter as explained 
in Ref. [SU] we get the ensemble-averaged G{V,T). 

The results for the Ohmic regime (Fig. [2]) were ob- 
tained by choosing a very large uj = 40 to ensure V ^ T. 
We analyzed two different system lengths: L — 30 — 7.5a 
and L = 50 = 12.5a. For L = 30 we generated an en- 
semble of 20, 000 samples and for L = 50 we used 10, 000 
samples in order to average out the mesoscopic fluctua- 
tions. (Actually, using 500 samples would give results of 
comparable quality.) Figure [2] clearly demonstrates more 
than a decade of the APLB of Eq. ^ near the inflection 
points of the curves. Note that this point is located at a 
lower temperature for the longer sample. 

Having determined the range of T where we get the 
Ohmic APLB, we proceeded to analyzing the non-Ohmic 
behavior of the system in this range of temperatures. To 
this end we fitted the results for higher V to Eq. Q. 
For L = 30, Fig. [l] we found a good collapse in both 
the Ohmic and non-Ohmic regimes. All curves in Fig. [l] 
were cut at V = 2, since at that point the curves were 
beginning to saturate as they approached the maximum 
current possible in the system. The collapse obtained for 
L = 50 (not shown) was equally good. The quality of 
our data collapse matches or exceeds that in the experi- 
ments [8l [T3HT9] . The values of the fitting parameters a, 
/?, and 7 have already been discussed (see more below). 

Let us now examine how the APLB we have found nu- 
merically can be understood in the light of established 
theory of ID VRH. According to this theory, transport is 
characterized by several regimes. At low T, the conduc- 
tance of the ensemble is dominated by rare paths with 
nearly equidistant sites, see Fig. [Sjb). This regime was 
studied in Refs. \28 \ \29 \ and [33l for two intermediate sites 
and in Refs. [26] and [27| for a chain of many sites [34] . 
Adopting the derivation in Ref. to the ID case we 
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FIG. 3. (a) A typical hopping path through the wire (thick 
line). The thin line represents electrochemical potential (b) 
A rare path [26H29| made of equal-length hops. Here 5e ~ NT 
and 5x ~ Na, where N is the number of hops. Weakly varying 
ri{x) corresponds to the Ohmic regime. 



can show that at a given T the main contribution to Gq 
comes from the chains of TV = \f2Lj\a hops, where A 
is the solution of the equation A ~ \w{\lLgT). In addi- 
tion, [MlISniEHl 



dlnGn 
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At T = Tinf, we have A ~ 1, which yields Eq. 

For T > Tinf, the system enters the regime where 
the transport is limited by rare highly resistive links — 
"breaks" — on the optimal path [5BH38j . As a result, 
the Ohmic conductance, which can be derived from the 
formulas of Ref. [351 obeys Eq. ([T]) with the coefficients 
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The concavity of the In Gn vs. In T curve is opposite in 
the two temperature ranges, which creates the inflection 
point, see Fig. [2] 

The non-Ohmic transport is also characterized by an 
^-shaped curve of In/ vs. InV, with its own inflection 
point. For example, at T 2> Tint, the theory ^39^ predicts 
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By the argument similar to that used in the Ohmic 
regime, (3 + 1 in Eq. ^ is determined by the maximum 
slope, i.e., the derivative of In / with respect to \nV: 



dlnV 



(15) 



We see that (3 + 1 (x T~^/^ is not a constant but de- 
creases with r, in qualitative agreement with experi- 
ments [3 dSl [17] . This explains why the data collapse 
onto the universal curve of Eq. Q can be achieved only 
in a limited range of T. 

At T = Tinf, Eq. (IT4| is at the border of its validity. 



Hence, Eq. ( 15 ) gives only the order of magnitude esti- 
mate, /3+1 < a, which is the first part of Eq. Finally, 
to get 7 we note that the crossover to the T-independent 
behavior in Eq. Q takes place at 7F ~ 27rT. On the 
other hand, according to Eq. ( 14 ), this occurs at uj ^ um 
where V/Tq ~ LT /oTq. Combining these expressions, we 
recover the second part of Eq. ^ . 

Formulas ([5])-([6]) predict numerical values and rela- 
tions among a, /?, and 7 that are in agreement with 
most of the cited experiments [HI ITBHTO] . Additionally, 
they provide a way of estimating the localization length 
a. For example, taking parameters a = 4.3, /3 = 2.1, 
7 = 0.25, T ~ l/xm of a representative MoSe2 nanowire 
from Ref. [ini we find a ^ 40 nm for this sample (W3). 
Other samples measured in that work showed a oc 1/\/M 
scaling with the number of transport channels M . In our 
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model the same scaling occurs if a cx M, as in a weakly 
disordered quasi-lD metal. In such a system a can be en- 
larged by applying an external magnetic field [40]. This 
is one convenient way to further test our model experi- 
mentally. Alternatively, it may be possible to vary the 
disorder strength and therefore a by electrostatic gating, 
while monitoring the predicted trends in a, /?, and 7. 

Another model we considered in search for the APLB 
was the interrupted-strand model (ISM) |4jjl42j. Therein 
a metallic wire is divided into segments by randomly po- 
sitioned impurities of tunneling transparency 1, 
which turn it into a chain of weakly coupled quantum 
dots. In the simulations we studied wires with Ni — 50 
impurities of strength s = 4. While we did observe the 
APLB in such wires {a = 3.75, /3 = 1.6, 7 = 0.15), the 
quality of the data collapse was not as good as in Fig. [T] 
Otherwise, the results were similar [15] . 

Note that the VRH in the ISM is analogous to multi- 
dot cotunneling in a granular metal. The latter also 
leads, in all spatial dimensions, to the power-law con- 
ductance behavior with a = f3 — 2Nd — 4, Nd being the 
number of dots involved in one cotunneling event |44| . 
Hence, the APLB is not uncommon in the VRH regime. 

The authors arc grateful to K. Matveev and M. Bock- 
rath for discussion of the results and to M. Di Ventra for 
providing access to computer resources on which a part of 
the calculations was carried out. This work is supported 
by the Grant NSF DMR-0706654. 
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